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ABSTRACT 

The problem of finding a function which, in addition to being zero 
outside a specified range in x-space, has its spectral power well con- 
fined to a certain range in k-space is solved numerically. Properties 
of the solutions are also discussed. 
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It is common practice to filter a set of data before performing a 
spectral analysis. Ibis is done either to remove unwanted noise or in 
hope of enhancing a weak signal thought to be present in the data. As 
both the filtering and spectral analysis are almost always performed on 
a digital computer, the filter "window" is effectively zero outside some 
range due to the finite precision of computer arithmetic. Therefore we 
will consider a filtering function g(x) which is defined to be zero 
except for x e[-L, L] . We have chosen to consider a symmetric range to 
simplify the analysis; however, we will see that the results are easily 
generalized. It is well known that if a function is zero outside [-L, L] 
there is no finite K for which its Fourier transform is zero outside 
[-K,K], Nonetheless, the aim of the filtering usually is to remove or 
retain some range of frequencies. Hence we want to consider a filter 
which is zero outside [-L, L] but also has its Fourier transform concen- 
trated in the range [-K,K]. Before going further we must choose some 
measure of how well a function's Fourier transform is concentrated in 
[-K, K], The measure we will use is the fraction of the spectral power 
contained in [-K,K]. For a given K and L we may then ask if there 
is a function which is optimal in the sense of being zero outside [-L, L] 
and having the largest possible fraction of its spectral power in [-K,K], 
We choose the following definition of the Fourier transform: 

( 1 ) 


_ 1 5 ? .4 W 

'g(k) = (2jt)" 2 f e" g(x)dx 


Then from the Fourier theorem we also know that 
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Then the quantity we wish to maximize is 



We will restrict our attention to real valued filter windows. The problem 
as presently posed is indeterminate without some normalization constraint, 
so we make the obvious choice 


00 



g 2 (x)dx = 1 
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We see now that our choice for a measure of compactness in k- space makes 

generalizing the range in both k- and x- space very easy. Translating 

the range in x-space AL say to [-L+ Ai*> L + AL] will multiply the 

Fourier transform by e lk ^ L , leaving the spectral power unchanged. 

x Akx 

Similarly multiplying the filter by e will translate the Fourier 

transform in k-space. Maximizing (3) is equivalent to finding the 
eigenfunction associated with the largest eigenvalue of the integral 
equation 


a f(x) = ~ C sin(x - x ) f ( x /) dx / for | x | < KL (5) 

where f(Kx) = g(x) and C£ is the fraction of the spectral power of g 
contained in [-K,K]. This result is derived in appendix I. In appendix II 
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it is shown that if f satisfies (5) then its Fourier transform satisfies 


a f (( kl) k) = i f sin(k - ^ f((KL)k') dk' . (6) 

J k — k 

-KL 

f cannot be its own Fourier transform since f(x) is defined to be zero 
outside [ -KL, KL] while (6) is valid for all k. Nevertheless if the 
largest eigenvalue of (5) is non- degenerate, f and its Fourier transform 
are simply related. 

We now make use of the fact that the kernel of (5) has a diagonal 
expansion in spherical Bessel functions (Jackson I962) : 


sin(x - x' ) 


x - 


= (2n+ l)j n (x)j n (x / ) 

n=0 


( 7 ) 


Now if we define « (x) = ( 2 n+ l)^\j (x) and expand the eigenfunctions 
T n n 

in these (non- orthogonal) basis functions we can rewrite (5) as 


.KL ec 


a f(x)=a ^ a n cp n (x) = | f ^ ^( x )cp n (x') ^ . (8) 

n =0 -KL n =0 m =0 

y*KL 

If we define the inner product (f,g) = I f(x)g(x)dx we may use the 

-KL 

independence of the spherical Bessel functions to write 


“ a n = \ £ ( V%. )a m • (9 

m=0 

Since our range of integration is symmetric, the inner products (cp n > 9 m ) 
are zero unless m+ n is even. Thus the problem separates naturally 
into finding the largest eigenvalue and associated eigenfunction for 
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odd and even functions separately. Therefore the remarks at the end of 
appendix II apply to all of the eigenfunctions of this kernel. We will 
leave implicit the fact that there are two separate problems to be 
solved. 

Now we may attempt to approximate the eigenvalues and eigenfunctions 
by truncating the spherical Bessel function expansion at some point. 

Thus we need a reasonably efficient method of calculating the spherical 
Bessel functions and the inner products of (9). The spherical Bessel 
functions are easily calculated by "backwards recursion" ( see, for 
example, Abramowitz and Stegun 1964a). This method is particularly 
suited to this case since the value of several spherical Bessel functions 
will be needed for the same argument. The inner product can be relatively 
easily calculated from a relation given by Abramowitz and Stegun (1964b) 
from which one can show 


J = iT+'m+T 



(x)j m (x) 


+ 


CO 


2 E 

k=l 


J n+k 



( 10 ) 


Therefore we need only evaluate the spherical Bessel functions at the 
end points of the interval of integration to find the required inner 
products. For this calculation standard library subroutines were used 
to evaluate the largest eigenvalue and associated eigenvector a^ of (9) 
truncated at N (for details of the algorithm used see Wilkinson and 
Reinsch 1971a, b,c). Then similar results were calculated for ( 9 ) 
truncated at N+ 1. When the difference between succesive approximations 
to the eigenvalue and the eigenfunction evaluated at zero and the end 
point of the range was sufficiently small, the approximations were 
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accepted. The fraction of spectral power in [-K,K] is plotted as a 
function of KL for both the even and odd cases in figure 1. The best 
eigenfunctions for several values of KL are plotted in figure 2 with 
the axes scaled so that the functions are non- zero on the same range 
and have the same value at zero. 

The solution of this problem may also be of interest to workers 
writing computer codes that use a k-space representation truncated to a 
finite number of harmonics to solve a differential equation in x- space 
( see, for example, Buneman 3-97^) • In this application a function in x-space 
is represented by a finite range in k-space (actually by a finite number 
of harmonics) but it is desired to restrict the function to a suitable 
range in x-space. Thus we have found a k-space form factor that is in 
some sense optimal for this problem. It is also possible that the solution 
to the similar problem in higher dimensions may be of interest. The 
diagonal kernels for the two and three dimensional cases are derived in 
appendix III. Ihe expansion of the eigenfunction in (8) also leads to an 
expansion of the eigenfunction in Legendre polynomials as shown in 
appendix IV. 
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APPENDIX I 


We wish to maximize 



(i.i) 


under the constraint 


f g^(x) dx = 1. 

“L /w 


Now we can rewrite (i.l) using 


the Fourier theorem and the requirement that g be real and zero 
outside [-L, L] 



K L L 

U 

•C-K •i-L -L 


g(x)g(x / ) 


ikx 


e 


-ikx 7 

e 


dx dx 7 


dk 


( 1 . 2 ) 


Interchanging the order of integration and performing the integral over 
k yields 


a = si - n - g( x) g( x 7 ) dx dx 7 

ki L 

r L 2 

We want to maximize (X under the constraint I g (x)dx = 1. 

•i-L 

we introduce a Lagrange multiplier and take the variation 


( 1 . 3 ) 


Therefore 



d.*) 
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For an extremal an arbitrary variation must be zero so we may set the 
expression in curly brackets equal to 0. 


1 r sin K(x-x') , , s . , , s „ I I ^ r 

— I - — - g(x')dx' = y g(x) for |x| < L 

J X - x / 

— ±J 

multiplying both sides by g(x) and integrating identifies 
and setting f.(kx) = g(x) results in equation 5 . 


( 1 . 5 ) 
as a 
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APPENDIX II 


We will demonstrate the validity of equation (6). First we define 
a function H(z) and write down the integral equation for its Fourier 
transform 


a H(k) 


1 

it 



3 in(k-k') g(k , )dk , 
k - k' 


(ii.i) 


Now we take the inverse Fourier transform of both sides which gives 


a h(z) = (2it) 


03 

J 


^ihz dk 



H(k')dk' e ikZ dk 

k - k' 


(II. 2) 


Interchanging the orders of integration on the right hand side and 
performing the integral over k yields 


a h(z) = ^fej^(cr( z+ l) " a(z-l)j f H(k') e lk z dk' ( IT . 3) 

-KL 

where u(z) is the Heaviside step function. Now we replace H(k' ) by 
its definition and interchange orders of integration obtaining 


a H(z) = i i|cr(z+ 1) - ff(z - l)j j* 5 H(z' )dz' f e -ik (z - z ) dfc / . (II.4) 

-KL 

Performing the integration over k' and noting that any function 
satisfying (ll,4) will be zero outside [-1,1] we may write 
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(ii. 5 ) 


a 




i 


z - z‘ 


Thus letting z = — we see that 
KL 


satisfies the same integral 
equation as f(x) so that f(k) satisfies equation (6). We can say 
more in the case that the eigenfunction is either even or odd. Since 
the square integrals of a function and its Fourier transform are equal 
we can in these cases write 


= (£ f e 19 f ext (( KL ) k ) 


(n.6) 


where f is the analytic continuation of the eigenfunction and 
0 = 0 or jt for even eigenfunctions, while @ = ± jt/2 for odd eigen- 
functions . 
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APPENDIX III 


In this appendix we briefly outline the method of solution for the 
two- and three-dimensional analogs of the problem treated in the text. 
The problem to be solved in two (three) dimensions is the following: to 
select from the set of functions which vanish outside a circle (sphere) 
of radius L that member which has the largest fraction of its spectral 
power contained in a circle (sphere) of radius K in k- space. From 
the symmetry of the problem we assume that the solution is cylindrically 
(spherically) symmetric, i.e. a function of the radius only. The two 
cases follow. 

A. Two Dimensions 

For a function f(r) vanishing for r > L, the two dimensional 
Fourier transform is a function of the magnitude of k only and is 
defined as 

f'(k) = ^ rdr f d<|) e _lkr COS< t ) f(r) = f rdr J Q (kr)f(r) , (ill.l) 

4 ^0 u 

where J (x) is the ordinary Bessel function of the first kind. The 
inverse transform is given by 


f(r) 



kdk J^(kr)f(k) 


From (ill.l) and (ill. 2) we see that 



(III. 2) 
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where 


p(r) = n/t"' f(r) and p(k) = Vk' f(k). We wish to find f (or f) 
such that the quantity 


a = ir dk r 2 (k) (III. 

is maximized. As in the main part of the text, we set N = 1 as our 
normalization constraint. From (lll.l), (ill. 3) then becomes 

(in. 4) 

where 


a = 


M 


dr, r(i , )H(r,r')r(r / ) 


H(r, r' ) 



kdk J Q (kr) J Q (kr' ) 


= K 




r' J^( Kr' 



(HI. 5) 


(in. 6) 


To maximize CL subject to the normalization constraint, we proceed as 
in appendix I to find that p(r) must satisfy the homogeneous integral 
equation 


«r(r) = j H(r,r')r(r')dr' . (ill. 7) 

4 ) 

The maximum value of 0! is then the largest eigenvalue of (ill. 7 ) 
and the desired function P(r) is the corresponding eigenfunction. To 
calculate CL and find the eigenfunction P as we did in the text 
requires that we find an expansion of the kernel H(r,r') analogous 
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to ( 7 )* For these purposes the expression of H(r,r') via (ill. 6) is 
not particularly useful so we return to (ill. 5 ) and note that 


1 /*2ir 

r o( kr ) J o( kr ') = 2S| J 0 ( kR ) d t 


(III. 8) 


2 2 2 

where R = r + r' - 2rr / cosf. But 


S' 


kdk J Q (kR) 


I J i (KR ) 


= 2 


y; (n+ 1 ) 

n=0 


Vl (Kr > J nt-l (Kr,) gln(nH)» 


r/ 


sinty 


(111,9) 


( III. 10) 


(cf. Gradshteya and Ryzkik 1965 a, b). From (ill .8)-(lII .10), (ill. 5 ) 
becomes 


H(r,r') 


CO 

2 ^(2irtl) 

n=0 


J 2nt l^ Kr ^ 




\lr^ 


(III. 11) 


Proceeding in a manner similar to that presented in the text, we then 
find from (ill. 7 ) and (ill. 11) that "r(r) can be written as 


r(r) = V?Z a ni (Kr) 


(ill. 12) 


n=0 


where 


<|> n (r) = V2(2n+1)' J2irf-l (lCr) 

and the column matrix a.^ satisfies the eigenvalue equation 
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(in. 13) 


aa = V' A (KL) a 
n / a nm m 

n=0 

wi th 


A 

nm 


(KL) 


KL KL 

= f -T = 2 V(2rH-l)(a^x5 r f . 

Vo •'o 


From here on a method analogous to that employed in the one-dimensional 

case can be used to evaluate the entries in the matrix A and to 

nm 

solve the truncated form of the matrix equation (ill. 13 ) for the largest 

eigenvalue and corresponding eigenvector a . From (ill. 12). then, 

n 

our desired function f(r) is given by 


f ( r) = A r(r) = 1 = f) 2 ^ a n 

n=0 n=0 




B . Three Dimensions 

Again in this case, the function f(r), vanishing for r > L, has 
a three-dimensional transform, which is a function of |k| (= k) only, 
given by 


/ 




d^ 


-ikr(_i 

e 



rdr sin kr f ( r) . 

( III .14) 


The inverse transform is 


rf(r) 



kdk sin kr f ( k) 


(III. 15) 
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From (iII.lA) and (IH.I 5 ) we see that 



( r)dr 



r 2 (k)dk s N 


where p(r) = rf(r) and T(k) = kf(k). 
The quantity to be maximized is 



(III. 16) 


Again setting N = 1 as the normalization constraint and using (III.14) 
in (ill. 16 ) we find 



r(r)H(r,r')r(r') 


where 


( hi. 17 ) 


H(r, 



dk sin kr sin kr' 


1 

jt 


sinK( r-r / ) 
r-r' 


1 

jt 


sinK( rt r ' ) 
r+r' 


( hi. 18) 


The extremal condition for OC with the normalization constraint is 
then the homogeneous integral equation 


ar(r) = f H(r,r')r(r' )<**“' , 

4 ) 


or, from (lll.l8). 


OT(r)=-jJ dr' slnK (^- r/ ) F ( r /) - i J dr' g^ r - r ') T(-r'). (m.19) 


r-r' 


U 



Thus if we define p(r) for r < 0 according to 

r(r) = -r(-r) , (III. 20) 

we see that p( r ) must satisfy 

ar(r) = i f L dr' sinK ( r - r ') r (r') . (ill. 21) 

51 L 

We note that (III.21) is exactly the same equation as that for the 
one- dimensional case ( 1 . 5 ). However we have the added restriction^ 

(ill. 20), that p(r) be odd. Thus the solution in this case is the 
same for the one-dimensional case, except that we must find the largest 
eigenvalue corresponding to an odd eigenfunction r(r). The solution 
for f(r) which has the largest fraction of its power spectrum confined 
to a sphere of radius K in k- space is then given by 

f(r) = 7 r(r) . 

The largest value of Oi as a function of KL for three dimensions is 
plotted in figure 1 as the curve corresponding to the best odd function. 
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APPENDIX IV 


We have from equation (6) the Fourier transform of an eigenfunction 


/. 


H 

03 

(Wl) a j 2n (KLk) j 

II 

'm 

V* 

1 

1 2 &n ' 

| A ( 


1 

[* 1 J 

n=0 , 

/ W 3) s j 2iH . 1 (KLk)j 


where the upper quantities in curly brackets refer to even eigenfunctions, 
the lower quantities to odd eigenfunctions. We have the general form 
for the Fourier transform of a spherical Bessel function from Abramowitz 
and Stegun (1964c) 


V k) “ (if ( ' l)n P n (k) (IV ’ 2) 

where P (k) is the Legendre polynomial of order n. When we take the 
inverse Fourier transform of both sides of (lV.1) we find 

*<*> ■ *(A) 4 1, <-«” 

(<‘«> 4 '„(il) 

We see that the value of the first neglected term gives an upper bound 
on the error in truncating (lV.3). 
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FIGURE CAPTIONS 


Figure 1. Fraction of spectral power confined (a) and log^(l-Ci^) as 
a function of KL. Curve 1 is the fraction of spectral 
power contained in [ ~K, K] for the best even eigenfunction on 
[-L,L] as a function of KL. Curve 2 is the same quantity 
for the best odd eigenfunction. Curve 3 is the logarithm of 
(1-a) , the fraction of the spectral power outside [-K, K], 
for the best even eigenfunction on [-L, L] as a function of KL. 
Curve 4 is the same quantity for the best odd eigenfunction. 

Figure 2. Eigenfunctions associated with the largest eigenvalue for 
various values of KL. The eigenfunctions associated with 
the largest eigenvalue plotted for KL varying from jr/2 
( uppermost, curve) to 4 jc (lowest curve) in increments of 
jt/2. As the functions are in these cases symmetric, they are 
plotted only for x > 0. The ordinate is arbitrary with all 
functions normalized to have the same value at x = 0 . 
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FIGURE 2 
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